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It is demonstrated that the TaN molecule is the best candidate to search for T,P-violating nu¬ 
clear magnetic quadrupole moment (MQM), it also looks promising to search for other T,P-odd 
effects. We report results of coupled-cluster calculations of T,P-odd effects in TaN produced by the 
Ta nucleus MQM, electron electric dipole moment (EDM), scalar—pseudoscalar nucleus—electron 
interactions, also of the molecule-axis hyperfine structure constant and dipole moment. Nuclear 
calculations of ^®^Ta MQM are performed to express the T,P-odd effect in terms of the strength 
constants of T,P-odd nuclear forces, proton and neutron EDM, QCD parameter 9 and quark chromo- 
EDM. 


INTRODUCTION 


During last years low-energy experiments on heavy 
atoms and diatomic molecules containing heavy atoms 
have proved to be very important for the search for 
New physics beyond the standard model [J-Q- The best 
limit on the electron electric dipole moment (eEDM) 
was obtained on the ThO molecular beam in [sj. The 
experiment was also sensitive to other effect that vi¬ 
olates time-reversal (T) and parity symmetries (P) - 
scalar—pseudoscalar nucleus—electron neutral current in¬ 
teractions 0. If a heavy atom has nuclear spin I> 1/2 
one can expect that the molecule will also be sensitive to 
other T,P-odd effect, the T,P-odd interaction of nuclear 
magnetic quadrupole moment (MQM) with electrons IB- 
Ill- Important ideas on the subject were suggested in 


Mrl 

papers |5l4ll|. 

There is a number of systems on which experiments to 
search for T,P-odd effects have already been conducted 
or suggested and which are inv estigated theoretically and 

YbF [i,|Ii®jThO 0, 
’ si, PbF [iilai, RaO 
4^ . etc.). 

Recently, TaN molecule was suggested as a new system 
to search for the T,P-odd MQM of the tantalum nucleus 
0 , where the molecule was marked as quite promising 
candidate to search for T,P violation in nuclear sector 
using molecules. Indeed, ^®^Ta nucleus is stable and its 
MQM is strongly enhanced due to the collective effect [ 1 ] , 
large nuclear charge Z of Ta leads to a higher than Z'^ 
electron enhancement of the T,P-odd effects iSiini, 
TaN in ^Ai state has D-doublets which allow full polar¬ 
ization in a small electric field and cancellation of many 
systematic errors due to opposite signs of the effects on 


exp erimentally (HfF+ 
[iaialll, ThF+ Ji3Jp, 
^^1, RaF PtH+ 



the doublet components [lll,|43, very small magnetic mo¬ 
ment due to cancellation of the orbital and spin contri¬ 
butions 0 and long lifetime of ^Ai electron state (43 |. 

TaN molecule has the triple bond order, i.e. higher 
than in some other recently experimentally considered 


systems to search for T,P-violation effects, e.g. ThO 
molecule 0. This suggests a rather complicated elec¬ 
tronic structure in terms of electron correlation effects. 
Up to now the molecule was rather Mor investigated 
both experimentally and theoretically M, ,0. There are 
no experimental data about the molecule-frame dipole 
moment, hyperfine structure, etc. Therefore, the aim of 
the present paper is to perform the first reliable ab initio 
treatment of the TaN electronic structure in ^Ai state 
including different T,P-violating effects. Also we pro¬ 
vide the values of the molecule-frame dipole moment and 
molecule-axis hyperfine structure constant which can be 
important for experimental planning and determination 
of the quantities and acceleration of the further experi¬ 
mental treatment of the system. To solve the problem we 
have developed pure ab initio scheme of the calculation 
where electron correlation effects were considered within 
the most accurate all-order coupled-cluster method, in 
which we took account of up to quadruple cluster ampli¬ 
tudes. 

Spin-dependent T,P-odd Hamitonian can be expressed 
in the following form 0 : 


TJT.P ^ ^SP ^MQM 


The T,P-odd scalar—pseudoscalar nucleus—electron 
Hamiltonian with a characteristic dimensionless 

constant /csp is given by [ 0 : 

= *^^A:sp7o75PAr(i’), (2) 

where Z is the heavy nucleus charge, Gf is the Fermi- 
coupling constant, 70 and 75 are the Dirac matrices and 
Pn(j) is the nuclear density normalized to unity. For the 
interaction of eEDM with inner molecular electric field, 
E, one has 


= 2de 


0 0 \ 

0 crE ) 


( 3 ) 
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where is the value of eEDM, cr are the Pauli matri¬ 
ces. Hamiltonian of interaction of MQM with electrons 
is given bv 47. 5^1 


^MQM ^ _ 


M 3 [a X r]^rk 

2I[2I-iy^2 r5 


(4) 


where Einstein’s summation convention is implied, a. are 
the 4x4 Dirac matrices, ^^,ris the displace¬ 

ment of the electron from the Ta nucleus, I is the nuclear 
spin, M is the nuclear MQM, 




3M 


2 /( 2 /- 1 ) 
2 


Tr 


i,k 


Ti^k — li^k + Ikh — + l) ■ 


(5) 

( 6 ) 


In the subspace of ±17 states (D = J • n|4'), J is the 
total electronic momentum, 4' is the electronic wave func¬ 
tion for the considered ^Ai state of TaN) the expressions 
(I2I3I4I) are reduced to the following effective molecular 
Hamiltonians, correspondingly 


HeS = W^SP kspS' ■ n. 

(7) 

= Wd d,S' ■ n, 

(8) 

o-MQM WmM 

- 2 /( 2 /- 1 )* ’ 

(9) 


where n is the unit vector along the molecular axis C 
directed from Ta to N, S' is the effective electron spin |5l| 
defined by the following equations: >= DID >, 

S'±|H = ±1 >=0 [13,141, 5'=|H|=1. 

To extract the fundamental parameters A:sp, de, M 
from an experiment one needs to know the factors Hsp, 
Wd and IEm, correspondingly, which are determined by 
the electronic structure of a studied molecular state on a 
given nucleus (discussed in Refs. [13,147, M, M|): 




fcsp 


( 10 ) 




( 11 ) 


= ( 12 ) 

i \ t / 4 

Note, that a parameter known as the effective electric 
field acting of unpaired electrons, /?eff = is often 

used. 

Eor a completely polarized molecule the energy shift 
due to eEDM, scalar—pseudoscalar nucleus—electron 
neutral current, and MQM interactions are 


(13) 


(5sp = kspWsp^, (14) 


dM{J,F) = {-iy+^C{J,F)MWMn , (15) 


C{J,F) 


J 2 J\ 

-non) 

2 ( I 2 I 

\-I 0 I 


(2J+1) 


J I F 
I J 2 


,(16) 


where (...) means elements with 3j—symbols and {...} 
with 6 j—symbols [s^, F is the total angular momen¬ 
tum and J is the number of rotational level. Note, that 
both Sd and (5sp are independent of J and F quantum 
numbers, whereas 6m depends on them. Also, in op¬ 
posite to i/gg and //gjj Hamiltonians //^^'^ has non¬ 
zero off-diagonal matrix elements on J quantum num¬ 
ber (between different rotational levels). This should 
be taken into account when mixing of different rota¬ 
tional levels become significant. In Eq. m this ef¬ 
fect is neglected. For ^®^TaN (1=7/2) and ground ro¬ 
tational level J=1 Eq. (flbll gives the MQM energy shifts, 
|(5(J,/’)|, equal to O.IOTIEmM, 0.1431TmM, 0.05WmM 
for F = 5/2, 7/2,9/2, correspondingly. 

To compute hyperfine structure constant A|| on ^®^Ta 
in the ^Ai electronic state of ^®^TaN molecule the fol¬ 
lowing matrix element can be evaluated: 

= i^)> (17) 

7 ' 7. / C 


where ^xa is the nuclear magnetic moment of a Ta isotope 
with spin /. 


NUCLEAR MAGNETIC QUADRUPOLE 
MOMENT 

The main contribution to MQM is produced by the 
nucleon-nucleon T,P-odd interaction which exceeds the 
nucleon EDM contribution by 1-2 orders of magnitude 
[^. In a spherical nucleus MQM is determined mainly 
by a valence nucleon which carries the nuclear angular 
momentum / Q: 

M = [d —2 ■ I0~‘^^ri(^ — q)(e ■ cm)]Ap(2/ — l)t/, (18) 

where tj = 1 for / = I + \/2 and tj = —1/(1 + 1) for 
I = Z — 1/2, I and I are the total and orbital angu¬ 
lar momenta of a valence nucleon, ij is the dimension¬ 
less strength constant of the T,P-odd nuclear potential 
r]G/(2^/^mp)(a ■ Vp) acting on the valence nucleon, p 
is the total nucleon number density, the nucleon mag¬ 
netic moments are Pp = 2.79 for valence proton and 
Pn = —1.91 for valence neutron, qp = 1 and = 0, 
Xp = h/rripC = 2.10 • 10“^^ cm, the contribution of the 
valence nucleon EDM d was calculated in Ref. @ . 

In deformed nuclei MQM has a collective nature and is 
enhanced by an order of magnitude [ 3 | (comparable to the 


6d = deWdn, 
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enhancement of an ordinary collective electric quadrupole 
moment). An estimate of the collective MQM can be 
done as follows. In a deformed nucleus the strong field 
splits orbitals with different absolute values \Iz\ of the 
projections of the angular momentum on the nuclear 
symmetry axis. Sum over all Iz gives zero contribution 
to MQM. However, due to the difference of energies for 
different \Iz\ some of the Iz orbitals are vacant, and the 
remaining contribution is not zero. As a result, the MQM 
of a deformed nucleus in the “frozen” frame (rotating to¬ 
gether with a nucleus) may be estimated using the fol¬ 
lowing formula Q: 

= E I.,l)nil, h, 1), (19) 


where the sum goes over occupied orbitals, 
Iz,l) is given by Eqs. ((TH|) and (O, 
Tzz = ‘2Iz — ^I{I + 1)? n{I,Iz,l) are the orbital 
occupation numbers, which may be found in Ref. 

The sum over a complete shell gives zero; therefore, for 
shells more than half-filled it is convenient to use hole 
numbers in place of particle numbers, using the relation 
M|“si®(hole) = (particle). 

The nucleus ^®^Ta has the following occupation num¬ 
bers: 10 neutron holes in orbitals [li,Iz] = [P3/2,±3/2], 
[zi 3 / 2 ,±13/2,±11/2], [h 9 / 2 ,± 9 / 2 ,± 7 / 2 ], arid 9 pro¬ 

ton holes [(^ 5 / 2 , ±5/2], [/iii/ 2 , ±11/2, ±9/2], [(i 3 / 2 , ±3/2], 

[57/2,7/2]. 

The MQM in the laboratory frame M = Miab can be 
expressed via MQM in the rotating frame ([T^ : 


^lab ^ 
(l.l77p 


(/ + l)(2/ + 3) "" 


— 0.9r/„) • 10 ^^(e-cm^) 


— {S.Odp + 2.3d„) • 10 ^^cm. 


( 20 ) 


where / = 7/2 is the nuclear spin of ^®^Ta. 

The T,P-odd nuclear forces are dominated by the ttq 
meson exchange Q. Therefore, we may express the 
strength constants via strong ttNN coupling constant 
g = 13.6 and T,P-odd ttNN coupling constants corre¬ 
sponding to the isospin channels T = 0,1,2: rjn = —gp = 
5 • 10®g(5i-l-0. 4 ^ 2 — 0 . 2 ^ 0 ) (see detailes in Hi). As a re¬ 
sult, we obtain 


M{ 9 ) = -[g{gi ± 0.452 - 0 . 250 ) • 8 • lO'^^e • cm^. (21) 


Possible CP-violation in the strong interaction sector is 
described by the CP violation parameter 0. According 
to Ref. ISa 550 = —0.375. This gives the following value 
of MQM for 181 Ta: 


M{0) = —5 • 10 ^85 . e . cm^. 


( 22 ) 


Almost the same final results for M (5) can be obtained by 
using recently calculated relations of 50 and 51 constants 
by using chiral perturbation theory |56l ] and lattice data 


for the strong part of the proton-neutron mass difference 


57| which give updated value of 550. 


Finally, we can express MQM in terms of the quark 
chromo-EDM and dd using the relations ggi = 
4..10i'i(J„ - Jd)/cm, 550 = 0.8 • 10^^{du ± d„)/cm [bi]: 


M{d) = —3 • 10 11 (d„ — dd) ■ e • cm. 


(23) 


The contributions of dp and dn to MQM in Eqs. (ED-El) 
are from one to two orders of magnitude smaller than the 
contributions of the nucleon CP-odd interactions. 


ELECTRONIC STRUCTURE CALCULATION 
DETAILS 


To obtain the electronic state-specific parameters lUsp, 
EeS, Wm and Ajj described by Eqs. (I11I10I12I17P we 
have performed a series of calculation with using the two- 
step procedure to study the relativistic four-component 
electronic structure in the vicinity of the Ta nucleus 
n in, . For this the space around the given heavy 
atom, Ta, is divided into the valence and core regions. At 
the first step inactive core electrons are excluded from 
molecular calculations using the generalized relativistic 
effective core potential (GRECP) method [ 6 ^| 6 lj. The 
approach allows one also to take into account the con¬ 
tribution from Breit interaction and finite nuclear size 
6 ^ 63 j. After this stage we obtain wavefunction that is 
very accurate in the valence region but has incorrect be¬ 
haviour in the core region. The correct four-component 
behaviour of the wave-function in the core region is re¬ 
stored at the second step using the procedure [32], [blj 
based on a proportionality of valence and virtual (unoc¬ 
cupied in the reference Slater determinant) spinors 
in the inner-core regions of heavy atom. The splitting of 
solution of the four-component relativistic full-electron 
problem on two consequent steps allows one to consider 
high-order correlation effects (see below) that are impor¬ 
tant for reliable and accurate calculation of the properties 
that cannot be measured experimentally. 

The applied computation scheme includes the follow¬ 
ing steps: (i) Consideration of leading correlation and 
relativistic effects within the relativistic two-component 
coupled-cluster with single, double and perturbative 
triple cluster amplitudes. The correlation calcula¬ 
tions included valence and outer-core electrons of Ta 
(bs^bp^bs^bd^) and N {\s^2s^2p^), i.e. 20 electrons. For 
Ta we used the basis set consisting of 15 s—, 10 5 —, 10 
d—, 5 /— and 2 5 — type uncontracted Gaussians [65|. For 
N we used the aug-ccpVQZ basis set with removed 
two g-type basis functions. Below this basis set will be 
called MBas. (ii) Consideration of correction on basis set 
enlargement. It was calculated as a difference between re¬ 
sults of 20 -electron scalar-relativistic calculation within 


the CCSD(T) method in enlarged basis [67| and basis 
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used at step (i). (iii) Consideration of high-order correla¬ 
tion effects. They were calculated as a difference between 
results of 20-electron two-component calculations within 
the couple cluster with single, double, triple and per¬ 
turbative quadruple cluster amplitudes, the CCSDT(Q), 
and the CCSD(T) calculation with using a compact basis 
set [Slii. ( iv) Consideration of contribution of sub¬ 
outer-core 4s^4p^4d^^4f^^ electrons of Ta. Their contri¬ 
bution was obtained as a difference between the results of 
52-electron and 20-electron 2-component CCSD(T) cal¬ 
culations [g^. Core electrons of Ta (Is — 4/) in steps (i), 
(ii) and (iii) were excluded from correlation treatment by 
the 60-electron GRECP. For step (iv) we constructed the 
GRECP version for 45 explicitly treated electrons of the 
Ta atom and a very small 28-electron inner core. 

The coupled-cluster calculations were performed us- 
ing the MRCC 0 [nl code interfaced to DIRAC code 
[72l| . Scalar-relativistic calculations were performed us¬ 
ing CFOUR program package [t^- Restoration of four- 
component electronic structure was performed using the 
code developed in [Hjlllilzll and interfaced to the above- 
mentioned packages. 

Equilibrium Ta—N distance and other spectroscopic 
parameters were obtained by approximation of poten¬ 
tial energy curve calculated within the 20-electron two- 
component CCSD(T) method in MBas basis set. 


RESULTS AND DISCUSSIONS 


smaller than Eeff in ThO [2^ [s^ HI L The value of 
the Wm parameter estimated in Ref. 4^ (« 1) agrees 
with the value obtained in this paper. Thus the calcu¬ 
lated large value of the Wm parameter conhrms that the 
^®^TaN molecule with stable ^®^Ta isotope having nucleus 
spin I> 1/2 probably the most promising candidate to- 
date among considered heavy-atom diatomic molecules to 
search for nuclear MQM. According to the results listed 
in table El and our previous error analysis from Ref. 
for a comparable situation (ThO in ® Ai state) we expect 
that the theoretical uncertainty of the calculated char¬ 
acteristics is less than 7%. Direct experimental check of 
Eeff, ITsp and Wm parameters is impossible, however, 
one can in principle seize here the opportunity of mea¬ 
surement of the hyperfine structure constant in ® Ai state 
of TaN which can be performed later. As it was argued 
earlier [^, [s^ the “equal-footing calculation” of the 

the hyperfine structure constants can provide a very im¬ 
portant (though indirect) test for the Eeff, Wsp and Wm 
parameters. The calculated value of the molecule-axis 
hyperfine structure constant is given in table El 

One can express the MQM energy shift, 
(—1 )^+'^C'(J, E)MWmO in terms of the fundamen¬ 
tal CP-violating physical quantities dp, dn, 0 and du,d 
using Eqs. 1)2012212311 . For the lowest rotational level, for 
which the coefficient |C'(J=1,E=7/2)| = 0.143 reaches 
a maximum value, we have 


0.143 WmM = - 


1025(4.5dp -f 3.5d„) 


e • cm 


^Hz (24) 


The equilibrium internuclear distance calculated for 
the ®Ai state of TaN is 3.19 a.u. which agrees well with 
the experimental datum see table ID In calculations 
of the parameters under consideration we have set R(Ta- 
N) to 3.20 a.u. 

TaN has triple bond order (one cr-bond and two tt 
bonds) provided mainly by d-electrons of Ta and p- 


0.143ITmM = -7.7 • 1O^°0 • pRz 


electrons of N [TGKTOj. Unpaired electrons in ®Ai state 
are non-bonding electrons localized on Ta. They deter¬ 
mine the configuration, where tr is mainly the 6s 
atomic orbital of Ta and 5 is mainly bd atomic orbital of 
Ta. 

The calculated values of Egff, IUspi Wm and A|| are 
given in table [TTl The effective electric field is rather 
strong (1.5 times more than that in HfF+ 14, [l^ and 
similar to other considered transition elements com¬ 
pounds Eeff(PbF)[3^ and WC [s^) but about two times 


0.143IUmM = -4.6- 


W^du - da) 


cm 


/iHz 


(25) 


(26) 


< 8.6 


The current limits on dp, |d| and |d„—dd| (Mj 
10“^®e-cm, |0| < 2.4 • 10“^°, \du—dd\ < 6 • 10“^”^ cm 
correspond to the shifts |0.143 WmM\ < 40 /iHz, 18 /rHz 
and 28 /rHz, respectively. Currently the best limit on 
the energy shift produced by the T,P-odd effects (elec¬ 
tron EDM and scalar—pseudoscalar nucleus—electron in¬ 
teractions) in ® Ai state of ^®^ThO is 700 /iHz Q and is 
expected to be improved by an order of magnitude over 
the next five years [a [13 due to experimental develop¬ 
ments. Such a progress with the ^®^TaN experiment on 
the ® Ai state can result in comparable improvement with 
limitations on the T,P-odd effects in the nuclear sector 
as well. 


TABLE I. Equilibrium internuclear distance Re, harmonic vi¬ 
brational wavenumber cUe and vibrational anharmonicity oje^e 
for the ® Ai state of TaN. 


Method 


Re, a.u. CUe, cm ^ UJeXe, Cm 


20e-2c-CCSD(T), this work 3.19 
Experiment, M 3.20 


1028 


3.5 
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TABLE 11. The calculated values of the molecule-frame dipole moment (d), effective electric field {Ees), parameter of the 
T,P-odd scalar—pseudoscalar nucleus—electron interaction (Wsp), parameter of T,P-odd MQM interaction {Wm), hyperfine 
structure constant (A||) of the ®Ai state of TaN using the coupled-cluster methods. 


Method 

d, 

Feff, 

IFsp, 

Wm 

All 


Debye 

GV/cm 

kHz 

10^^ Hz 

MHz 

20e-2c-CCSD 

4.70 

37.5 

34 

1.15 

-3082 

20e-2c-CCSD(T) 

4.81 

35.6 

32 

1.10 

-3015 

correlation correction 

-0.04 

0.2 

0 

-0.01 

-32 

basis set correction 

-0.03 

-0.4 

0 

0.00 

-33 

core correction 

0.00 

-0.5 

-1 

-0.01 

-51 

FINAL 

4.74 

34.9 

31 

1.08 

-3132 


“ Magnetic moment of Ta was set to 2.371 um [80I. I81ll . 
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